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Abstract 

A lattice approach is used to describe the mechanical interaction of a corroding rein- 
forcement bar, the surrounding concrete and the interface between steel reinforcement 
and concrete. The cross-section of the ribbed reinforcement bar is taken to be circular, 
assuming that the interaction of the ribs of the deformed reinforcement bar and the sur- 
rounding concrete can be captured by a cap-plasticity interface model. The expansive 
corrosion process is represented by an Eigenstrain in the lattice elements forming the 
interface between concrete and reinforcement. Several pull-out tests with varying degree 
of corrosion are analysed. The numerical results are compared with experiments reported 
in the literature. The influence of the properties of concrete are studied. The proposed 
lattice approach offers insight into corrosion induced cracking and its influence on bond 
strength. 



1 Introduction 



Corrosion of reinforcement involves the transformation of steel into rust, which is an 
expansive process [S]. If the expansion is restrained, it induces radial pressure in the 
confining material. For reinforced concrete, the radial pressure and accompanying trans- 
verse tensile stress may cause cracking [5]. Cracking is not desirable because it reduces 
the anchorage capacity of the reinforcement P, [191 E] ■ Most of the anchorage capacity of 
deformed reinforcement is provided by ribs on the surface of the bar, which resist the slip 
between concrete and reinforcement by transferring inclined radial forces into the con- 
crete [2^. The capacity of the concrete to resist these forces can be significantly reduced 
by corrosion-induced cracking. Consequently, there is a considerable interest in develop- 
ing models of corrosion-induced cracking, which can quantify the influence of cracking on 
the bond capacity of reinforced concrete. 

The mechanics of corrosion-induced cracking and its influence on bond properties are 
complex. For example, at the microscale, the compaction and penetration of rust into 
pores and micro-cracks takes place [2X1127] . Modelling at this scale is challenging, since the 
microstructure of concrete is complex. Therefore, most research is concentrated either at 
the mesoscale, where the interaction of the ribs of the reinforcement and the heterogeneous 
concrete are described explicitly, or the macro-scale, where these interactions are described 
by an interface constitutive model. At the macro-scale, the formation of rust is modelled 
by an expansion of the interface, which results in macroscopic cracking and subsequent 
spalling. Many of the models at the macro-scale proposed in the literature include the 
effect of corrosion-induced cracking on bond by reducing the bond strength of the interface 
between concrete and steel [19] . In these models, the relationship between the amount of 
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Figure 1: Three-dimensional modelling of reinforced concrete: Concrete, reinforcement 
and bond between concrete and reinforcement are considered as individual phases. 

rust and the reduction of bond strength is determined empirically. Thus, these models 
are of limited validity for the prediction of the influence of corrosion on bond. Other 
models describe the expansion of the rust, the radial pressures and the transverse stresses 
on the concrete explicitly [T7] . These models have the potential to establish an analytical 
relationship between the expansion of rust, cracking and spalling. They can be combined 
with realistic bond models [20] , so that the influence of corrosion-induced cracking on the 
bond capacity can be predicted. However, this modelling framework is computationally 
demanding, since it requires three-dimensional modelling of the mechanical response of 
the concrete, the bond between reinforcement bar and concrete, and the reinforcement 
bar itself, as shown in Fig. [ij 

In view of these difficulties, discrete methods appear to offer a favourable alternative ap- 
proach, as they are known to be efficient for the modelling of displacement discontinuities 
at interfaces [18]. Discrete methods can be subdivided into two categories: particle models 
and lattice models. In particle models, the arrangement of particles can evolve, so that 
former neighbours may become separated. Therefore, such models are suitable for de- 
scribing processes involving large deformations. On the other hand, in lattice models the 
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connectivity between nodes is unchanged during analysis, so that contact determination 
is not required. Consequently, lattice models are mainly suitable for analyses involving 
small strains [151 [Ml II]- Their potential to model corrosion-induced cracking and its 
influence on bond is assessed in the present study. 



2 Lattice modelling approach 

Lattice approaches have been used successfully in the past to model the failure of concrete 
[21113]. In these approaches, the lattice elements do not represent the meso-structure of the 
material, but instead discretise the continuum. For example, Bolander and his coworkers 
have accurately reproduced analytical solutions for elasticity and potential flow problems 
[2H1 E] ■ Lattice models can also incorporate constitutive models, formulated in terms of 
tractions and displacement jumps, as commonly used in interface approaches for concrete 
fracture p]. These models yield element size- independent descriptions of crack-openings 
[12]. The spatial arrangement of the lattice elements and their cross-sectional properties 
is based on Delaunay and Voronoi tesselations of a set of random points placed in the 
domain [2H] • The random placement of nodes reduces the influence of the discretisation on 
the fracture patterns, as observed for other fracture approaches [131 IB] • For concrete, we 
adopt a damage-plasticity constitutive model following the work presented in [HI [T31 [10] . 
For the bond between reinforcement and concrete, a non-associated plasticity interface 
model is proposed, which is conceptually similar to the model developed by Lundgren 
[20] . This model consists of a Mohr-Coulomb friction law combined with a compression 
cap. For the steel phase, an elastic constitutive model is used. 

For each two-noded lattice element, a local co-ordinate system is introduced (Fig. [2]): the 
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axis n is aligned with the element axis and axes s and t are aligned with the two principal 
directions of the cross-section of the lattice element. Each node has six degrees of freedom 
(three translations and three rotations) which are used to determine the displacement 
jump at the centroid c of the element cross-section in the local coordinate system by rigid 



body motions as 



Ur = BUr. 



where Uc = {ucn,Ucs,Uct,(pn,4's,0t}'^ are the displacement and rotation discontinuities at 
the point c and Ue = {ui,Vi,Wi, 0i, z/'i, ^i, ^2, "^2, '?^2, 02, "^2, ^2}'^ are the degrees of freedom 
at the two nodes. Furthermore, the matrix B in Eq. Q is 
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In Eq. ([2]), Cs and Ct are the eccentricities between the midpoint of the lattice element 
and the centre c in the directions s and t, respectively. In addition, h is the length of the 
element and A is the cross-sectional area. Furthermore, /p is the polar second moment of 
area and Ji and I2 are the two principal second moments of area of the cross-section. In 
the local coordinate system, the stiffness matrix is: 
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Figure 2: 3D lattice element. 

where D is the 6x6 constitutive matrix of the material. In the following sections, the con- 
stitutive models for bond, concrete and reinforcement are described for the displacement 
discontinuities Uc = {ucn,Ucs,Uct}'^ ■ 

2.1 Constitutive model for the bond between concrete and re- 
inforcement 

The interface response between reinforcement and concrete is characterised by displace- 
ment jumps which are related to tractions via an interface constitutive model. The 
three-dimensional displacement jump Uc = {wcn, Mcs, Wct}"^ is transformed into strains 
£ = {^nj^sj^t}^ by means of the interface thickness h as 



The thickness of the interface h is chosen to be equal to the length of the lattice element, 
which for the bond model crosses the interface between the reinforcement steel and the 
concrete. The strain s is related to the stress cr = (o"n, cxs, o"t)^ by the elasto-plastic 
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. transition 




Figure 3: Yield surface: Mohr-Coulomb friction law combined with a cap. 



stress-strain relationship 



Cr = De — — Scor) 



(5) 



where De is the elastic stiffness, £p = (^pn, ^ps, ^pt)"^ is the plastic strain and e, 



EcoT, 0, 0) is the Eigenstrain used to describe the expansion associated with the corrosion 



process (see Sec. 2.2). The elastic bond stiffness matrix is 



^ ^ 

7b^b 



7b^b 



(6) 



where i?b is the normal stiffness and 7b is the ratio of shear and normal stiffnesses. The 
yield surface of the plasticity model consists of a Mohr-Coulomb friction law combined 
with an elliptical cap. The shape of the cap surface is adjusted so that a smooth transition 
between the two surfaces is obtained (Fig. |3]). This combination was initially proposed 
by [25] for a circular cap and further developed by [7]. The yield function / depends on 
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the normal stress cXn and the shear stress norm = A/cr|~T~of 



as 



CTq + ObCTn 



(cTn + /c - af 



PI 



PI 



if (Tn < (TnO 



(7) 



where ctb is the friction angle and /c is the compressive strength of concrete. Furthermore, 



/3bttb/c 



ab/3b + + 



where /3b is the ratio of the short and long axes of the cap elhpse (Fig. [3]). At the point 
where the two parts of the yield surface meet, the normal stress is 



/^bftb^l + /^b«b 



(9) 



The rate of the plastic strains in Eq. ([s]) is 



OCT 



(10) 



where g is the plastic potential and A is the plastic multiplier. In the present study, g is 
chosen to be very similar to the yield function /. The only difference is that Ob is replaced 
by a flow inclination ipT^ so that the magnitude of the normal plasticity strain generated 
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during shear loading can be controlled. Thus, the plastic potential is 



9 = < 
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/3bV^bv/l + /3bV 



(13) 



The plasticity model is completed by the loading and unloading conditions: 



/ < 0, A > 0, A/ = 



(14) 



This plasticity bond model is similar to the one developed by Lundgren [20] for bond 
between concrete and reinforcement, but there are several differences. Here, a smooth 
transition between the cap and the frictional law is introduced, so that a special vertex 
stress return algorithm in the transition region is obviated. Also, the response is perfectly- 
plastic. Thus, the calibration is different from the one used for the model proposed by 



Lundgren. 



9 



2.2 Model for corrosion between concrete and reinforcement 



The effect of corrosion is modelled by an Eigenstrain ^cor in Eq- which is determined 
from the free expansion of the corrosion product Ucot as 



(15) 



-cor , 

h 



where h is the length of the element. For the bond model, this is the length of the element 
across the concrete-steel interface. This approach has been shown to give results which 
are independent of the element length |9]. The free expansion Ucor is determined from the 
corrosion penetration Xcor, which is the thickness of the layer of steel that is lost during 
the corrosion process (Fig. lik). The percentage of steel loss p is 



— ^ X 100 = — ^ — ' X 100 16 



where is the cross-sectional area of the uncorroded reinforcement bar, AV"s is the cross- 
sectional area that is lost during the corrosion process and is the diameter of the 



reinforcement bar. Solving Eq. (16) for a;cor gives 



Xr 



0/2 (l - - P/lOo) (17) 



The transformation of the steel layer Xcor into rust is assumed to be accompanied by a 
cross-sectional area expansion 

AKor = A,n,AK (18) 
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Figure 4: Representation of the corrosion process as an expansive layer of rust: a) steel 
loss, b) unrestrained expansion of rust 



assuming an unrestrained expansion (Fig. |4p). Here Acor is a model parameter. The 
cross-sectional area of the rust can be written as 



AKor = TT (0 + 2neor)' /4 - TT (0 - 2x,,,f /4 



(19) 



where Ucov is the expansion of the unrestrained corrosion layer as shown in Fig. Hb. Equat- 



ing Eqs. (18) and (19), substituting for AV^ from Eq. (16) and solving for Ucor gives 



Mcor = a/0V4 + {(pXcor " xl^^) (Acor - 1) - 0/2 



(20) 



The bond model, which considers corrosion, has seven parameters, namely i?b, 7b, /c, 
«b, Phi i'h and Acor- The calibration of these parameters is discussed later. 



2.3 Constitutive model for concrete 

The constitutive model for concrete is based on a damage-plasticity framework. The 
strains, which are determined from the displacement jumps as discussed earlier, are related 



11 



rp 

to the nominal stress cr = ((Xn, Us, (Xt) as 



<T = (1 -a;)Dc(s-£p) = (1 - (21) 

where w is the damage variable, De is the elastic stiffness, £p = (^pn, ^ps, £^pt) is the 
plastic strain and a is the effective stress. The elastic stiffness has the same format as 
the one for the bond model. The two parameters denoted and 7c control the Young's 
modulus and Poisson's ratio of the material. The plasticity model for concrete is based 
on the effective stress and is thus independent of damage. 

Again, the plasticity part is described by the yield function (Eq. [t]), the flow rule (Eq. 11 ) 



and the loading unloading conditions (Eq. 14 ) as described for the bond plasticity model. 



The yield function, shown in Fig. |5]for / = 0, which is a function of two stress variables 
cTn and d"q = a/o"s + cT^ , is defined as 



f={ 



2-2 , o"c(/c-«c/3c/t)_ , _2 2a2/Jt + (1 - ac/3c) ^ fc - aMt 

at.at. + 2 (Jn + if (Tn > — 7^ — 

l + ac/9c ^ l + ajc I + 01 J, 

I g /c - «c/^c/t _ -2 (1 - ac/^c) /c - 2Qc^c/c/t - /c - Qc/^c/t 

I (51 ^ /32 (1 + + ^ + (1 + ^^^j 1 ^ ^ ^^^^ 

(22) 

where ft is the tensile strength, /c is the compressive strength, and ac and /3c are the 
inclinations shown in Fig. |5} The plastic potential g is almost identical to the yield 
surface, except that the inclination ac is replaced by ipc- 



The damage parameter in Eq. (21 ) is determined by means of the damage history variable 



= {Spn) (23) 
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where (.) denotes the McAuley brackets (positive part of operator). The relation of this 
history variable to the damage parameter is derived from the response in pure tension for 
which (Tn > and Uq = 0. For this stress state, the crack opening is defined as 



Wr 



/i(£:pn + a;(£:n-^pn)) 



(24) 



where h is the length of the lattice element. The damage law is based on an exponential 
stress-crack opening relationship in the post-peak regime for the normal stress component. 
It has the form 



/texp 



Wf 



(25) 



where Wf is a model parameter, which is related to the local fracture energy Gp of the 
concrete lattice elements as Wf = Gp/ft- The normal stress component is also prescribed 



in Eq. (21) as 



(Tn = (1 - OJ) Ec {Sn - ^pn) 



(26) 



In Eq. (26), the elastic strain — Spn in the post-peak regime can be expressed as ft/Ec 



since the plasticity model is perfect plastic. Equating Eqs. (25) and (26), and using 



Eq. (23), a nonlinear equation for the damage parameter u is obtained as 



[1 — Lo) = exp 



h (/td + ft/Ec) 



(27) 



From this equation, the damage parameter u is determined iteratively by means of a 
Newton iteration. The eight parameters for the constitutive model of concrete are Ec, 7c, 
ft, fc, ac, /3c, ipc and Wf. 
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Figure 5: Yield surface: The yield surface is controlled by the tensile strength /t, the 
compressive strength /c and two parameters and /3c- 

2.4 Constitutive model for the reinforcement 

The material response of the reinforcement is modeled by the linear elastic stress-strain 
relationship 



where the two parameters of the elastic stiffness are denoted as and 7s. 

3 Comparison with experimental data 

The lattice approach described in the previous section is used to model the experiments 
reported in [19]. The geometry and loading setup of the experiments are shown in Fig. [6] 
Reinforcement bars (jzc = 13 mm) embedded in concrete cubes were initially subjected to 
corrosion and subsequently pulled out. In the analyses, it was assumed that the cubes were 
restrained on one face by friction-less supports, which were modelled by roller boundary 
conditions in the numerical analyses. All other faces are assumed to be traction- free. 
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Figure 6: Geometry and loading set-up for the corrosion pull-out tests reported by [T^. 
The reinforcement bar (diameter = 13 mm) is placed eccentrically. 

The concrete used in the experiments is characterised by a Young's modulus of 22.6 GPa, 
Poisson's ratio of 0.17, a tensile strength of 2.7 MPa and a compressive strength of 
24.7 MPa. The Young's modulus of the reinforcement is 183 GPa. The response of the 
concrete, reinforcement and bond between concrete and reinforcement is modelled using 
the lattice shown in Fig. [7| which consists of 9165 lattice elements. For the reinforcement 
and the interface between the reinforcement and the concrete the mesh is structured, 
but for the concrete itself the lattice is randomly orientated. The rectangles in Fig. [7] 
are the cross-sections of the lattice elements which represent the interface between the 
reinforcement bar and the concrete. For these elements, the element axis is normal to the 
interface. Therefore, the Eigenstrain normal component ^cor represents the strain in the 
radial direction. Furthermore, the cap-plasticity bond model for these elements is used 
to idealise the interaction of the ribs of the deformed bar and the surrounding concrete, 
which are not modelled explicitly. 

Three tests were analysed. In the first test, the reinforcement was pulled out without 
initial corrosion. In the other two tests, corrosion levels of p = 3.2 and 16.8 % were 
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Figure 7: A part of the mesh for the lattice analysis. 

reached before the puUout. In all three analyses, the load F was controlled by imposing 
an end slip, which is defined as the horizontal displacement of node A as shown in Fig. |6^. 

The proposed model for corrosion-induced cracking explicitly represents the three phases, 
(a) concrete, (b) reinforcement and (c) bond between concrete and reinforcement. Three 
constitutive models are required for this purpose, and inevitably this demands specifi- 
cation of several parameters; 17 in total. To determine these parameters, the following 
calibration strategy was adopted. Firstly, the two parameters for the reinforcement were 
calibrated based on a cube specimen so that the stiffness of the reinforcement was matched. 
This gives = 345 GPa and 7s = 0.065. The same specimen was used to calibrate the 
elastic parameters of the concrete model, which were determined to be Ec = 36.6 GPa 
and 7c = 0.175. In the next step, a uniaxial tension and compression test was carried 
out to determine the remaining parameters of the concrete model. Two uniaxial tests 
are not sufficient to determine 6 model parameters uniquely. Therefore, the parameters 
Oc = 0.5, f3c = 0.5 and ipc = 0.25 are assumed to be default values. The tensile and 
compressive strength were set to ft = 2.2 MPa and fc = 40 MPa. Furthermore, the 
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Figure 8: Comparison of predicted average bond stress-slip curves and experimental data 
reported by |l9j for three corrosion percentages p = 0, 3.2 and 16.8 %. 

parameter Wf = 0.045 mm was determined based on the assumption of a local fracture 
energy of Gp = 100 N/m. 

For the bond model, the elastic parameters were chosen as ii^b = '^EgEc/ {Ec + Es) and 
7b = 0.175. Here, is the harmonic mean of Eg and i?b assuming that half of the element 
crossing the interface is steel and the other half is concrete. The other parameters of the 
bond model were calibrated in two steps. First, the frictional parameter = 0.24 and 
dilation parameter ^/'b = 0.05 were calibrated so that the peak value of the experimental 
stress displacement curve for p = was obtained. The corrosion related parameter Acor = 
1.67 was determined by matching the peak stress of the experiment with 3.2 % of corrosion. 

The results of the analyses are compared to the experimental results in the form of average 
bond stress-slip curves shown in Fig. [8j Here, the average bond stress was determined as 
r = F/{n0i), where i = 60 is the embedded length (Fig. |6^). 
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The pre-peak regime of the load-shp curves obtained in the analyses is in very good 
agreement with the experiments for and 3.2 %, which is expected as these experimental 
results were used to determine the input parameters. However, the peak of the analysis 
for 16.8 % strongly underestimates the experimental result. Also, for all three corrosion 
percentages the post-peak response of the analyses is more brittle than the one observed 
in the experiments. The main failure mechanism observed in the analysis is the occurrence 
of splitting cracks which reduce the capacity of the specimen to resist transverse tensile 
stresses generated by the slip between reinforcement and concrete. This agrees with the 
observations made in the experimental study reported in [19]. For the analyses without 
corrosion, the crack patterns for slips of 0.38, 0.42 and 0.5 mm (marked in Fig.jsj) are shown 
in Fig. |9} Crack patterns are visualised as those middle cross-sections of lattice elements, 
in which the crack opening increases and exceeds 50 {im at this stage of the analyses. 
These cracks are called active. Just after the peak of the average bond stress-slip curve. 




Figure 9: Crack patterns for the pullout analysis for the corrosion-free case at (a) 0.38, 
(b) 0.42 and (c) 0.5 mm of slip. 

the concrete cover cracks at its thinnest section (Fig. |9^). With further slip, additional 
cracks initiate from the reinforcement and propagate through the concrete as shown in 
Fig. |9)d and c. The results of this pull-out test are further investigated by studying the 
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distribution of the shear stress along the interface. In Fig. lOi-c, the distributions for 



bond shps of 0.38, 0.42 and 0.5 mm are presented. Here, the coordinate x indicates the 
position along the reinforcement bar and a is the angle in the polar coordinate system 
defining the position at the interface between reinforcement bar and the concrete in the 
y-z plane (see Fig. |6b). It can be seen that, as the cracks propagate through the concrete. 
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Figure 10: Shear stress distribution along the interface for the three stages marked in 
Fig.i 



the shear stress becomes increasingly localized in the "crack-free" zones (at a = 0° and 
180°). 



For the corroded bars, the crack patterns are shown in Fig. 11 , at the end of the corrosion 



process before pull out. Again, only active cracks are shown. For both these cases, crack- 
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Figure 11: Crack patterns for the analyses with (a) 3.2 % and (b) 16.8 % corrosion before 
pullout. 
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Figure 12: Influence of fracture energy on the nominal bond strength. 

ing of the concrete cover occurs before pullout, which corresponds to the experimental 
observations reported in the literature 



Since the failure mode is dominated by the occurrence of splitting cracks, it is expected 
that a change in the capacity of the concrete to carry tensile stresses will strongly influence 
the observed bond behaviour. Therefore, a parametric study was carried out for all three 
cases, assuming three different fracture energies of the concrete lattice elements, namely 
Gp = 100, 200 and 400 N/m. The peak bond stresses obtained from these nine analyses 



are presented in Fig. 12 It can be seen that the fracture energy strongly influences the 



bond strength, particularly for highly corroded bars. For instance, an increase of the 
fracture energy from 100 to 200 N/m results in an 4.75 fold increase in nominal bond 
strength for a corrosion percentage of 16.8 %. This sensitivity of the results with respect 
to capacity of the concrete to transmit transverse tensile stresses may well explain the 
observed discrepancy between the numerical and experimental results, shown in Fig. [8} 
for the specimen with p = 16.8 %. For instance, in the analysis, the boundary at the 
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face of the specimen, at which the force was apphed, was assumed to be supported only 
in the direction of the apphed force. However, it is hkely that in the experiments some 
friction between the loading plate and the concrete specimen was present, which would 
have provided lateral resistance. Another simplification is the use of a constant expansion 
parameter Acor = 1-67 for varying steel loss. Experimental results reported in [27] show 
that rust products penetrate into the cracked concrete. This penetration is expected 
to be more significant in highly corroded bars. Furthermore, for highly corroded bars, 
compaction of the rust product as reported in [21] might play a role as well. 



4 Conclusions 

A lattice approach is used to describe the mechanical interaction of a corroding reinforce- 
ment bar, the surrounding concrete and the interface between steel reinforcement and 
concrete. The cross-section of the ribbed reinforcement bar is taken to be circular, assum- 
ing that the interaction of the ribs of the deformed reinforcement bar and the surrounding 
concrete is included in a cap-plasticity interface model. This lattice approach is capable 
of representing many of the important characteristics of corrosion-induced cracking and 
its influence on bond. The idealisation of the corrosion expansion as an Eigenstrain al- 
lows for the modelling of corrosion-induced cracking. The frictional bond law can model 
the decrease of the bond strength if the concrete is pre-cracked. Good agreement with 
experimental results in the pre-peak regime of the bond stress-slip curves was obtained 
but more studies are required to investigate the post-peak response of the bond stress-slip 
curves. It should be emphasised that the parameter controlling the amount of expansive 
Eigenstrain is not a material parameter of the rust, but a model parameter which takes 
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into account the effect of many micro and meso-structure effects, such as compaction of 
the rust and penetration of rust into cracks. In future studies, the dependence of this 
parameter on the amount of corrosion will be investigated. 
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